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We study the effect of crystal symmetry and step-edge diffusion on the surface current governing 
the evolution of a growing crystal surface. We find there are two possible contributions to anisotropic 
currents, which both lead to the destabilization of the flat surface: terrace current jt, which is parallel 
to the slope m = \7z{x, t) , and step current js , which has components parallel (jj ) and perpendicular 
(is") to the slope. On a high-symmetry surface, terrace and step currents are generically singular 
at zero slope, and this does not allow to perform the standard linear stability analysis. As far 
as a one-dimensional profile is considered, j^ is irrelevant and js suggests that mound sides align 
along [110] and [110] axes. On a vicinal surface, js destabilizes against step bunching; its effect 
against step meandering depends on the step orientation, in agreement with the recent findings by 
O. Pierre-Louis et al. [Phys. Rev. Lett. 82, 3661 (1999)]. 



I. INTRODUCTION 



The kinetic stability of a crystal growing homoepitaxially by Molecular Beam Epitaxy is determined primarily by 
the possible existence of a slope-dependent mass current j(m) along the surface, i.e. by a current which does not 
vanish in the limiting case of a constant slope m (m = Vz, where z{'r,t) is the local height) [Q-^. Such a current is 
generally ascribed to the so-called Ehrlich-Schwoebel (ES) effect at step-edges, which hinders interlayer diffusion |Q. 

On singular surfaces, experimental results (mainly on metal growth) show that the instability leads to mound 
formation and often to a coarsening process, where the typical size L of the mounds increases in time (generally 
with a power law: L{t) ~ i") |3|. The template of the mound structure is already formed in the early stages of 
growth (the so-called 'linear regime'), and here crystal structure should determine shape and orientation of mounds. 
For example, cubic crystals are characterized by a four-fold and a six-fold symmetry, respectively on (100) and 
(111) faces: experimental analysis by Scanning Tunneling Microscopy has indeed shown square based mounds on 
Fe and Cu(lOO) [^-Q] and triangular based ones on Rh and Pt(lll) ||,^. The relevance of the in-plane symmetry 
for the later stages of the growth process has been definitely proven by Siegert pO| , who has shown -through a 
continuum description of the surface- that unstable currents with different in-plane symmetries may give rise to 
different coarsening exponents n. For vicinal surfaces ES barriers at steps are known to stabilize against step bunching 
and to destabilize against step meandering ||ll|-]l^ . 

It is therefore extremely important to determine what are the microscopic mechanisms giving rise to slope-dependent 
currents j, what is the expression of j, and how lattice symmetry enters in it. One of the main results of the present 
paper is the finding of two contributions to the slope-dependent current, one due to terrace diffusion (jt) and one due 
to step diffusion (js)- Both contributions are singular a.t zero slope|^. This is at odds with the usual phenomenological 
expressions of jt, used in the continuum description of surface growth [|l0|,^6||l^] , which all reduce to the simple 
isotropic form jt — am in the small slope regime (to — > 0). We will sec that our expressions for jt and js remain 
anisotropic even in this limit, and that implies a singular behaviour in to = 0. Other important results concern the step 
current js, which is found to destabilize layer-by-layer growth against mound formation on a high symmetry surface, 
and step-flow against step bunching on a vicinal surface. Step-flow is stable or unstable against step meandering, 
depending on the step orientation. 

The destabilizing effect of ja on a singular surface has been observed independently by O. Pierre-Louis et al. jlj] 
and by Ramana Murty and Cooper jl^. The former have also studied analytically the effect on step meandering. 
Here we provide a unified treatment of these diverse effects within a continuum description of the surface, we predict 
the new phenomenon of step bunching induced by step currents, and we analyze the different anisotropic behaviours 
of and jt- 
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II. TERRACE CURRENT 



In Fig. |l| we draw a piece of a vicinal surface corresponding to a constant slope m, and a piece of a step. Once 
adatoms have landed on the surface, they perform a diffusion process until they stick to the upper or lower step. 
The attachement rate D' from below is considered extremely fast [D' / D — oo, D being the diffusion constant on the 
terrace); the rate D" from above defines the ES length £es = [D/D" — 1) (in units of the lattice spacing) |^,|l8|,|l9| . 
This should be compared to the diffusion length representing the minimal distance between nucleation centers on 
a high-symmetry surface Under the usual conditions of crystal growth we have £d ^ 1, while both the cases 

^ES *C (weak ES effect) and i?ES ^ (strong ES effect) may take place 



y 




FIG. 1. Sequence of equispaced steps, typical of a vicinal surface (left) and just a single step (right), taking into account the 
discrete character of the lattice. All the symbols are explained in the text. 



In one dimension (Id) we write the ES current (due to terrace diffusion) as jt = •mf{m?), and at small slopes (in 
the sense that m <^ 1/^d) we have the linear behaviour jt = am = /(0)m, with a — -F^es^d/[2(^es + ^d)], F 
being the intensity of the external flux (i.e. the number of particles landing on the surface per unit time and lattice 
site). In two dimensions (2d), if we neglect in-planc anisotropy we can generalize and write jt = rnf{rn?). Let us now 
discuss the microscopic origin of anisotropy and how it modifies jt. Throughout we will consider a (lOO)-surface with 
fourfold symmetry, and take the x {y) axis along the [100] ([DID]) orientation, denoting by x and y the corresponding 
in-plane unit vectors. The extension to other crystal symmetries is straightforward in principle. 

In the absence of surface reconstructions, terrace diffusion by itself is an isotropic process, at least in its continuum 
description. In contrast, the sticking of an adatom to a step depends on the microscopic environment, which depends on 
the step orientation. So, the ES barrier seen by an adatom approaching a step depends on the orientation of the surface 
and this dependence translates into an orientation-dependent ES length i?Eg = €es(^)i where = arctan(ma;/|mj,|) is 
the angle of the step relative to the x-axis. Assuming straight steps, the expression for a one-dimensional surface can 
be taken over, and for small slopes (m^D 1) we obtain 

^-^"^'^^"^ 2(^Es(^)+M "^- 

The coefficient a becomes independent of 9 only in the regime of strong ES barriers, ^es(^) ^ (in this limit, 
a — F£'^/2). For weak barriers in-plane anisotropy is therefore present even in the 'linear' regime m^D ^ 1 (a = 
^'^es(^)^d/2). Through the dependence of 9 on m, Eq.(|l|) is manifestly non-analytic at m = 0. 



III. STEP CURRENT 



Next show that crystal symmetry manifests itself also (and perhaps mainly) through step diffusion. Once adatoms 
have reached a step, they can diffuse along it at a rate Dg and stick to a kink edge (see Fig. |l|). Similarly to terrace 
diffusion, only if there is an asymmetry between D'^ and D'^ , a net step current js exists; the strength of the asymmetry 
determines an ES length along the step which will be called £k, the subscript k standing for kink. Step diffusion biased 
by kink barriers is similar to terrace diffusion along a one-dimensional surface |T^, but some differences are worth to 
be stressed. 

i) All the possible in-plane orientations 9 of the step, with the correct symmetries, should be taken into account, 
because -especially for a high-symmetry orientation- all the 9 are found on the same surface. This may be true also 
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for a vicinal surface, if steps are subject to a strong meandering ||13|. In particular, orientations corresponding to 
6 = and 9 = tt/A will be seen to behave in a qualitatively different way. ii) Adatoms arrive at the step at a rate 
Fs depending on the terrace size i. For equally spaced steps Fs = F£. However, since (in 2d) the surface current is 
defined as the number of atoms crossing per unit time a segment of unit length, orthogonal to the current, the actual 
expression for the current is obtained by multiplying the 'single-step current' by the number of steps per unit length, 
i.e. l/£. This factor cancels the factor i appearing in Fg, since the current is proportional to Fg as well, iii) A step is 
a one-dimensional object, and therefore it has a larger roughness than a two-dimensional surface. In the expression 
for the unstable current, the diffusion length gives the minimal distance between steps (in 2d) or kinks (in Id) along 
a high-symmetry orientation. In 2d, steps are created by nucleation and growth, and £d is generally given by an 
expression as lu « (D/F)'^ with the exponent 7 depending on the details of the nucleation process M. In Id, the 
corresponding expression £jj^ = (Dg/ Fs)^= should be compared to the distance between thermally excited kinks, and 
the smaller one (called id) be chosen. In most of our discussion we will assume that £d is sufficiently large so that 
double or multiple kinks can be neglected. 

The high symmetry in-plane orientations [100] and [110] for a step are fairly different in the mechanisms giving rise 
to a step-edge current. Along a [100] segment, step diffusion takes place between nearest neighbours lattice sites at a 
rate Dg^ and the analogy with a one dimensional surface is appropriate. In particular, an asymmetry in the sticking 
rates to a kink determines a net current along the straight segments of the step, i.e. in the x direction; when 6^0 
this current does not vanish and it has a component along the slope to, which will be seen to destabilize the flat 
surface. 

Conversely, along a [110] orientation, diffusion is a two-steps process and it is very much slower, because it requires 
detachment from a high coordination site. As a first approximation, it may even be reasonable to assume that no 
detachment at all takes place. This does not prevent a nonzero step-edge current, for the following reason. Along 
the [100] orientation, kinks are due to nucleation, or they must be thermally activated, because a kink increases the 
total length of the step. Along and close to the [110]_^orientation, the total length of the step does not depend on the 
specific sequence of [100] and [010] terraces (see Fig. ||), and therefore the step is rough even at zero temperature. The 
path for going from the origin O to P is equivalent to a directed random walk, where the asymmetry p between the 
probabilities p_ = (1 — p)/2 and p+ = (1 -|-p)/2 to go respectively in the x and y directions, is nothing but the tangent 
of the angle /3 = 7r/4 — 6 formed by the average orientation of the step with the [110] direction. Since step diffusion 
does take place along [100] and [010] segments, each step segment longer than one lattice constant contributes to the 
X and y components of the step current. This implies that jg is nonzero also for 9 = 7r/4: in this case jg is exactly 
parallel to m and it has a destabilizing character. In the following we will consider separately the cases of small 9 and 
9 close to 7r/4, and afterwards we will write down a general expression for jg, valid for any value of the angle 9. 

[010] [110] 




[100] 



FIG. 2. Step profile, when its average orientation is close to [110]. 

For the moment, we will suppose that the slope m = |to| of the surface is larger than 1/^d, i-e. we are in the 
'vicinal' regime. The step current jg can always be written as the sum of j| and j^, where j| = (jg ■ rn)rn/m'^ and 
j^ = (jg ■ rh±)rh±/m'^. The vector rn± = {—my,mx) is orthogonal to to. 

If we are close to the [100] orientation, the current jg is simply j[ioo] = jid(jns)x, where jidijng) is the usual 
unstable current for a one-dimensional surface whose slope is mg = tan0 = mx/\my\. For small 9, jid = dgvag, with 
ttg = [Fgj €)l].l?J\l{lk -I- ld)\- By decomposing j[ioo] along fh and toj^, we obtain 

-t jxdimg) ^ ^ 1 . . 

J[ioo] = 2 — V^^™' ~ mymj_\. (2) 
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The uphill component is (j[ioo] ' 'fn/m) « asrn^/rn^ > 0, for small TOs. The positive value of this component explains 
why step-edge current is enough to destabilize a flat, high symmetry surface. More details are given in Sec. 0. 

If we are close to the [110] orientation, as explained above, the current originates from entropic fluctuations around 
the straight step, which create segments along [100] and [010] directions. Each segment of length s contributes a local 
current proportional to (s — 1). Therefore, the components jx and jy of the step current along the x and y axes are 



simply proportional to the probabilities (per step site) 
horizontal and vertical direction, respectively: 



and pi_ to have a couple of consecutive step sites in the 



J[iio] = -f{p-x -ply) = F 



1-p 



F 



By using the relations p — tan{9 — tt/A) and iris = tan6', after some easy algebra we obtain 



J[110] 



F 



1 — TO, 



1 + |tos|to (1 + |tos|)^ m 



(3) 



(4) 



The expressions (H) and (^ are valid close to the [100] and [110] orientations: they can be generalized to any value 
of TOg, i.e. to any angle 6, by writing 



js = A{0)m/m + B{9)m^/m 



(5) 



Both A and B are periodic in 9, with period n/2 (see Fig. Their behaviours close to 6* = and 9 = n/A are 
derivable from j[ioo] (Eq. (^)) and j[iio] (Eq. (|4|)). More precisely, A{9) is always nonnegative, it has a minumum for 
9 = and a maximum for 9 = tt/A. The function B{9) vanishes at = 0, tt/A, it has a positive slope in = and 
a negative slope in 9 = n / A. These properties are all that we need in the following, and they do not depend on the 
specific model assumed to derive Eqs. (^,^), because they are mainly due to symmetry considerations. For example, if 
multiple kinks are allowed when the step is close to the [100] orientation, A{9) is no more zero at = 0, but A{0) > 
and 9 = is still a minimum. Finally note that in contrast to the terrace current jt |^], the step contribution (^) is 
independent^ of the surface slope m, i.e. the step distance, in the vicinal regime (to^d S> 1). 




n/4 Tz/2 ^0 7t/4\ /n/2 9 

FIG. 3. Plots of A{e) and B{e), which are defined in Eq. (i) 



N+-N. 
N 



Js,t 



Before going on, let us generalize the expression of js to any value of the surface slope m. In the limit to ^ 1/^d, 
both Js and jt go to zero, because contributions coming from steps and terraces of opposite sign compensate pp^. In 

a region of small slope, js.t = 

small slope 

(for js) or terraces (for jt) and N = 



, where N± is the number of positive and negative steps 

large slope 

iV_. Since (A^+ — N-)/N = to^d, in the small slope regime: jl = A{9)£j:,'rn 



and js = B{9)£]:,m±. For a generic slope, we can write 

j! = A{9)g{m'^)m 



jj- = B{9)g{m^)m^ 



(6) 



where g{m^) iu for to£d ^ 1 and g{m?) 1/to for ml-D 3> 1. The simplest function interpolating between the 
two limiting values is g{vn?) = Id/ ^Jl + ra^i'^, but its actual form does not need to be specified. 



■^Strictly speaking, a slope-dependence is maintained through the ^-dependence in Id,, but when i is small, (.d, should be 
replaced by the distance between thermally activated kinks (see the main text). 
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IV. STABILITY OF VICINAL SURFACES 



Let us now perform a linear stability analysis of a growing vicinal surface of average slope toq . The local height is 
z{x,t) = rriQ ■ X + e{x,t) and the local slope is to = toq + Ve. The evolution, as determined by the step current, is 
given by the equation dtz — —V • js- By using the general properties given above, for A and B, and working in the 
special cases = and 6*0 — tt/A, we obtain: 

dtz = -V • J. = -g{ml)[A{0o) + B'{0o)]dle - A{eo){d/dmo)[mog{ml)]dle (7) 

where d± and 9|| are directional derivatives perpendicular and parallel to to (i.e. parallel and perpendicular to the 
step). Thus the coefficient of d\ gives informations on step meandering, and that of on step bunching [ pl| . 

Since A{0) > and {d / dmo)[mog{mQ)] > 0, the current js (more precisely j|) has always a destabilizing character 
against step bunching; if multikinks are not allowed along the [100] orientation, ^(0) = and the effect is absent 
along the x,y axes. Also, as it may be expected, has no effect on this instability. 

Concerning step meandering, we must distinguish between 9o — and 9q = 7r/4, because B'{6o) has opposite signs 
in the two cases. For 9o = the derivative is positive and therefore steps along the [100] orientation are destabilized 
by step meandering. On the contrary, the evaluation of [A(7r/4) + i3'(7r/4)] gives, using Eq. (||), a negative result 
(= —F/^/2). Therefore, steps along the [110] orientations are stabilized against step meandering. 

Our conclusion regarding the [100] steps agrees with the analysis of Pierre-Louis et al. [^, while in the case of [110] 
steps they find stability (instability) for large (small) values of the kink ES length Ik- Our expression (^ is valid if 
adatoms are not allowed to turn around corners and this effectively sets ik — oo. Its generalization to any ik indeed 
shows that A{tt/A) becomes a minimum (and the quantity [A(7r/4) + i?'(7r/4)] changes sign) for £k < 2^*, where £* is 
the typical distance between corners along a [110] step. In the 'random walk' model for the step, £* — 2 and therefore 
the orientation [110] is stabilized unless £k is extremely weak. 

Ramana Murty and Cooper [^ have performed Monte Carlo simulations of a vicinal surface, with steps along the 
[100] axis. Step meandering is indeed observed, even if the terrace current jt is absent. Conversely, no step bunching 
seems to occur, suggesting that their simulations correspond to A{0) = 0. 

V. STABILITY OF SINGULAR SURFACES 

The analysis of a high-symmetry surface is complicated by the non-analytic behavior of js and jt in m = 0. In the 
small slope regime (m <C l/iu), js and j'^ become (see Eq. (|^)) 

jl = eDA{9)m jj- - £uB{9)mi_ . (8) 

It should be stressed that the singularity is physically justified, as we now try to argue. Close to an extremum of the 
profile, z{x, y) = zq + {ci/2)x^ + (c2/2)y^ -I- csxy, rh = {ax + c^y, C2j/ -|- cacc) and rus = (cir -I- c^)/\c2 + car], where 
r = x/y. The prefactors A and B are therefore manifestly non-analytic functions at x = y = 0. The reason is that 
close to an extremum, steps are closed lines and as the top (or the bottom) of the profile is approached, the step 
orientation is no more defined. The angular dependence of A and B also implies that the evolution equation does not 
become linear in the small-slope regime, and hence arbitrary profiles cannot be treated as superpositions of harmonic 
ones. 

The problem of non-analyticity does not appear when we consider a one-dimensional profile, i.e. a profile varying 
only in one direction (for example, z = z{x,t)), because the prefactors A and B are now constants^. This implies 
that the divergence of the current is easily evaluated: 

dtz = -V • {]s + ]t) ^ -V • [£YiA{9)rn + eoB{9)mj_ + a{9)m] 

= -[£uA{9o) + a{9o)]V^z = -v{9o)V''z . (9) 

The component of the step current parallel to the step (j^) does not contribute, because V • to_l = 0, while the 
component parallel to the slope {jl) destabilizes the fiat surface, analogously to jt- 



^The angle 6 may indeed take the values Oq and {6q -t-vr), corresponding to tan^ — ±|ms|, but because of the 7r/2 periodicity, 
^(^o) = A{9o + tt). This is no more true for the (111) surface of a cubic crystal. 
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The instability gives rise to pyramid-shape mounds, whose orientations 9i should be the most unstable ones, i.e. 
correspond to the maxima of v{9). In this respect, the step current favours the orientations forming 45° with the 
crystallographic axes, while the anisotropy induced by the terrace current depends on the microscopic details of the 
sticking processes. Since the presence of more kinks along the step should favour the descent of the adatom, it is 
reasonable to think that £es{(^) is maximum in 9 = 0,7r/2 |22|, and therefore the two contributions to i^(6') should 
compete. Close to 9 — 0, we have 



iy(9) = asTUs— + a{9) = a^r + a{9) 
m 



By using expression (Q) for a{9) and the expression of Os given above Eq. (|^), we obtain 



v'\0) = F 



2{i- 



ES 



(10) 



(11) 



Mounds will align along the crystallographic axes if ^{O) has a maximum in = 0, i.e. if i'"{Q) < 0. It is apparent 
that for a sufficiently large ES effect at steps this condition is not fulfilled, because the anisotropic character of jt 
disappears. On the other hand, £es should also not be too small, otherwise jt itself would be negligible. Therefore 
1/(0) will have maxima in 6* = 0, 7r/2 only if several conditions are simultaneously satisfied: i) £'-^^{0) must be negative^; 
ii) -^Es/^D should be neither too large nor too small; iii) £k/£d should be small. 

Refs. ||l5[ and ||2^ recently reported simulations on high-symmetry surfaces, taking into account step-edge diffusion. 
In both cases, if kink barriers are present the mound sides align along [110] and equivalent axes. Since in [ p^ there are 
no ES barriers at steps and in [p5[ the barriers are infinite and therefore isotropic, mound orientation is determined 
only by js, in agreement with our picture. 



VI. DISCUSSION 



Some aspects of the subject considered in the present paper have not been sufficiently clarified and they deserve 
further analysis. First of all, the non-analytic behaviour of the surface current remains problematic, because it 
implies that the continuum evolution equation 9(2: -|- V • j = is not defined at m = 0. As we have argued above, 
this non-analyticity is an inescapable consequence of the persistence of crystal anisotropy in the 'linear' regime of the 
instability; if it could be removed, e.g. through a more careful treatment of the interpolation between vicinal and 
singular surfaces, this would also imply that mounds are initially isotropic and develop their anisotropic shapes only 
in the nonlinear regime. It is however also conceivable that, as in the case of equilibrium surface relaxation below 
the roughening temperature p^ , the appearance of a singularity in the continuum evolution equation carries a real 
physical message: That the surface is not well described by a smooth function z{r,t) near its maxima and minima. 

A second important aspect concerns vicinal surfaces. We have seen that for [110] steps js has a stabilizing character 
with respect to step meandering and a destabilizing character with respect to step bunching. It would be interesting 
to evaluate quantitatively these effects and to compare them with the opposing effects of the terrace current. This 
comparison has been done for step meandering |lj] , and it seems that the effect of jg may dominate jt ■ At any rate 
the predicted step bunching instability should be clearly visible in simulations of models which have no ES barriers 
but only asymmetric sticking at kinks |l5|] . 

Finally, in this work we have not addressed the effects of crystal anisotropy on the smoothening terms in the 
continuum evolution equation, which are crucial in determining the actual length scale of the instability [^|T] . Under 
far from equilibrium conditions, the dominant smoothening mechanism is believed to be due to random nucleation [ p^ , 
which is manifestly isotropic; the anisotropy of the equilibrium step free energies will however be felt if detachment 
from steps becomes significant. 



*In simulations of solid-on-solid models the step edge barriers are often implemented such as to reduce the probability for 
adatoms to approach steps, rather than to descend from them |^,^. In this case the barrier at a [110] step may in fact 
(slightly) exceed that of the close packed [100] step [^ . 
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VII. CONCLUSIONS 



In this paper we have studied the different contributions to the surface current on a (lOO)-surface, which depend 
only on the slope m. Such contributions come from biased surface diffusion, both on terraces (jt) and along steps 
(is) J where the bias mechanism is an Ehrlich-Schwoebel barrier at steps and kinks, respectively. 

The expressions of js,t are relevant in two respects: they determine the linear stability of the flat surface and ~ in 
the case of unstable growth - they also determine the shape and the orientation of the emerging structure. 

The terrace current is parallel to the slope, while the step current has components parallel and perpendicular to 
the slope, because step diffusion takes place along the [100]- and [010]-segments that constitute the step. 

A first important result is that the anisotropic character of js and jt persists in the small slope regime rh 0: this 
means that they are both non-analytic at m = and consequently this does not allow a complete description of the 
evolution of a high-symmetry surface. 

Concerning the stability of a singular surface, js is found to be destabilizing because it has a positive component 
in the direction of rh. The most unstable orientations form angles of 45° with the crystallographic axes, while jt is 
usually (but not always, see footnote^) thought to select the x,y axes. If a competition exists js should prevail (see 
discussion below Eq. (|ll|)). 

The stability of a vicinal surface is more complex. The terrace current jt is known to stabilize against step bunching 
and to destabilize towards step meandering, whatever is the orientation of the surface. Surprisingly, the step current 
is instead found to generally favor step bunching (along the crystallographic axes js has no effect if multiple kinks are 
not present). Finally the effect of js on step meandering strongly depends on the surface orientation and the strength 
of the ES effect at kinks: [100] steps are destabilized, while [110] steps may be stabilized if the kink ES-barrier is not 
too weak. 
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support by the Alexander von Humboldt foundation. The work of J.K. was supported by DFG within SFB 237. 



APPENDIX A: LIST OF SYMBOLS 



z{x 
rfi = 
rfi I : 



t) 
Vz 

= {-my,m^) 

arctan ( ) 

- \m\ = \m^\ 
1/m 



mo 
Jt 

Js 

il' 

F 

Fs = 
D 

Ds 
D' 
D" 

D's 

^ES 

4 



Ft 



local heigth of the surface at point x and time t 
local slope 

vector orthogonal to the slope 

angle between the average orientation of a step and the x axis 

modulus of the local slope 
terrace size 

average slope of a vicinal surface 
terrace current 
step current 

step current parallel to the slope m 
step current perpendicular to the slope m 
intensity of the external flux of particles 
rate of particles arriving to the step 
diffusion constant on the terrace 
diffusion constant along a step 

attachement rate of an adatom on a terrace to the ascending step 
attachement rate of an adatom on a terrace to the descending step 
attachement rate of an adatom on a step to the ascending kink 
attachement rate of an adatom on a step to the descending kink 
ES length on a terrace 
ES length along a step 
diffusion length on a terrace 
diffusion length along a step 

the minimum between Iy)^ and the distance between thermally 
activated kinks 
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p = p+-p_ 



In the directed random walk picture of a step close to the [110] orientation, 
probabilities that the step goes in the x and y directions 
asymmetry in the probabilities p± 
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